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We derive finite-temperature sum rules for excesses in internal energy and in (volume-integrated) 
pressure arising due to presence of heavy quarks in SU{N) gluon plasma. In the limit of zero 
temperature our formulae reduce to the Michael-Rothe sum rules. The excesses in energy and 
pressure of the gluon plasma are related to expectation values of certain gluon condensates, and, 
simultaneously, to the heavy quark potential. The sum rules lead to a known relation between 
the internal energy and the potential, and to a new expression for the excess in the pressure. The 
pressure appears in the free energy as a generalized force associated with variations of the spatial size 
of the heavy-quark system. We find that the excess in gluonic pressure around a heavy quarkonium 
is always negative. Finally, we derive an exact equation of state that provides a relationship between 
the gluonic energy and pressure of heavy quarks. 
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Introduction. The most anticipated result of exper- 
iments on heavy-ion collisions is a possible creation of 
a deconfined state of matter in a form of quark-gluon 
plasma. Long time ago it was suggested that an enhanced 
dissociation of heavy quark-antiquark bound states could 
be a good signal of the deconfinement phase [ij. Theo- 
retical efforts aimed to study properties of heavy quarko- 
nia require a nonperturbative input in a form of quark- 
antiquark potentials at finite temperature 0, S] ■ A pow- 
erful numerical tool to study such potentials is based on 
first-principle calculations of Polyakov loop correlators in 
lattice simulations of QCD Q. The renormalized corre- 
lator of Polyakov loops provides us with the free energy, 
and - via thermodynamic relations - with internal energy 
and entropy of the heavy quarks Since the quarks 
are heavy particles they cannot be excited thermally, so 
that the term "free energy of heavy quarks" is usually 
understood as an excess in free energy of thermal gluons 
that appears due to presence of the strong chromoelectric 
fields of the heavy quarks [1] . 

It is clear that presence of heavy quarks affects ther- 
modynamics of the quark-gluon plasma. Despite the fact 
that the thermodynamic effect of one quark-antiquark 
pair is not an extensive quantity, a multiple production of 
quarkonia may provide a contribution to the bulk proper- 
ties of the plasma. In addition to the free energy, internal 
energy and entropy of the heavy quarks, we find that the 
quarks are also able to make contribution to the pressure 
of the system. The effect of pressure is not seen in a stan- 
dard approach because the pressure is usually associated 
with variations of the volume of the system, which does 
not enter the excess in the free energy of a finite-sized 
quark pair. However, below we show that the variations 
in size of a quark system do couple to the pressure. 

In the heavy-quark thermodynamics, the spatial extent 
of the quark pair can be treated as a external variable 
while the pressure enters a quantity that plays a role of 



the corresponding generalized force. From this point of 
view the renormalized quark potential can be associated 
with an excess in a generalized Helmholtz free energy. 

We use a sum rule approach that is generally known as 
powerful analytical for investigation of various nonper- 
turbative properties of QCD physics For example, 
in absence of the external sources the susceptibility of 
the trace of the energy-momentum tensor can be related 
via certain sum rules to the bulk viscosity in the quark- 
gluon plasma The susceptibilities of components of 
the energy-momentum tensor were shown to be related 
to (derivatives of) thermodynamic potentials 

We derive exact finite-temperature sum rules for ex- 
cesses in internal energy and in (volume- integrated) pres- 
sure in the presence of external heavy quarks. In the 
limit of zero temperature certain combinations of these 
sum rules reduce to the well-known action and energy 
sum rules derived by Michael and Rothe in a lattice for- 
mulation of the theory. The Michael-Rothe equations 
associate the quark-antiquark potential to an excess in 
chromoelectric and chromomagnetic condensates. His- 
torically, the first attempt to derive such sum rules was 
done in [Sj, but the derivation turns out to be correct only 
partly [9j. The rules were subsequently corrected [10| 
and extended [lH , and an important role of the confor- 
mal anomaly was stressed 1 121 . A lattice check of the 
sum rules was done in Ref. [13|. Below we use the new 
finite-temperature sum rules to study the heavy-quark 
thermodynamics further. 

Thermodynamics of SU(N) Yang-Mills theory. 

It is instructive to discuss first the thermodynamics of 
Yang-Mills theory without external quark sources. De- 
spite the derivation is well known, we repeat it for the 
sake of further convenience. 

The equation of state of the SU{N) Yang-Mills theory 
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can be determined from the energy momentum tensor 



= 2Tr 



(1) 



where G ^^^=0"^^^ =d fj,Ay - di,A^ + ig[A^, Ay] is the field 
strength tensor of the gauge field A^ = f^A"^, and t° 
(a = 1, . . . , N"^ — 1) are the generators of the gauge group. 

At a classical level the energy-momentum tensor ([1]) 
is traceless because the bare Yang-Mills theory is a con- 
formal theory. At the quantum level, however, the con- 
formal invariance is broken due to a dimensional trans- 
mutation. As a consequence, the trace of the energy- 
momentum tensor exhibits a conformal anomaly, 

%,ix) = ^G^,.(x)G;,(x) ^ ^dx) , (2) 
^5 9 

where C is the Euclidean Lagrangian of the theory: 



(3) 



Here and denote, respectively, chromoelectric and 
chromomagnetic field strength tensors squared. 
The Gell-Mann-Low /3-function in Eq. ^ is 



Pig) = JT] — = -5 ^0 
dm II 



(4) 



where fj, is the renormalization scale at which the Yang- 
Mills coupling constant g is defined. One- and two-loop 
perturbative coefficients of the /3-function ^ are 
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In the thermodynamic limit the energy of the system 
is E — eV where e is the energy density and V is the 
volume of the system (at the end of the calculation we 
set V^oo). The energy E and the pressure P are de- 
termined via the derivatives of the partition function 
Z — Z[V, T, gifJ.)] with respect to the temperature T and 
volume, respectively. 



E = T 



,91nZ 
ainT ' 



PV = T 



dhiZ 
d\nV 



(6) 



The logarithm of the partition function is a dimension- 
less quantity proportional to the free energy F = fV, 



InZ 



Thus, one gets on dimensional grounds 
/ d d d 



dhvV dlnT dlnfi 



InZ = 0. 



(7) 



(8) 



The free energy is a physical quantity, therefore its total 
derivative with respect to the scale /i should vanish: 



d 



dlnZ 



d In /i \ 9 In /i 



/3(g)^^lnZ = 0. (9) 



Here we used Eq. ^ and followed partly Ref. The 
identity d\nZ/dg = {V/T){2C/g) and Eqs. © , '© , ©, 
(El), ® give us the conformal anomaly relation 



Q = E-iPV= d 



y /?(g) 



E'^{x)+B^{x) 



(10) 



where we implemented the subtraction of the T — 
contribution, {{■ ■ ■))t — {■ ■ ■}t — {■ ■ ■)t=o UM- Equa- 
tion (fTO)) is the spatial integral of the conformally- 
anomalous trace ^ of the energy- momentum tensor ([l} . 

Equation (jlOp is implemented in most first-principle 
determinations of the gluonic equation of state from lat 
tice simulations of Yang-Mills theory (see, e.g., Refs. 
[l6jV Indeed, combining Eq. (O with © one gets: 
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E = F~T 
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df 



P 



F 
V 



(11) 



so that 6 
becomes 



d P[T) 



The left hand side of ^ 
while the right hand side can be de- 



ar 

termined numerically. Integration of this equation gives 
us full information on the thermodynamics of the theory. 

Thus, in thermodynamic limit of a uniform system - 
such as the vacuum of SU{N) gauge theory - the energy 
density and pressure are extensive variables that are re- 
lated to each other by thermodynamical identities. Equa- 
tion (|10p provides us with an additional relation between 
the energy and pressure of the system so that we have the 
complete system of two equations for two unknown quan- 
tities. However, the contributions to energy density and 
pressure coming from the external sources (for example, 
from heavy quark- ant iquark systems) cannot, in general, 
be described by extensive variables. Thus, the standard 
relation between energy and pressure should not work in 
these systems by default, and we need to determine these 
quantities separately. 

Fortunately, Eq. ([T0|) can be extended further by in- 
troducing separate renormalization scales fit and /x^ in, 
respectively, temporal and spatial directions. The proce- 
dure is well-known in the lattice gauge theory [l3, [3| ; 
and we reformulate it here in continuum terms for the 
future use. The analogues of Eq. ([8|) are: 



dlnZ 9 In Z 



dlnV dlnfi. 



-0, 



ainZ dlnZ 



dlnT dlnfit 



= 0, (12) 



while the requirement of the scale independence provides 
us with following counterparts of Eq. (|9]): 



dlnZ dlnZ 



d In fia d In fia 



E 

b=s,t 



dgb d\nZ 

d In fjLa dgb 



0. 



(13) 



Here a — s,t, and we always imply silently that at the 
end of all differentiations we set fit = fJ-s = The spatial 
and temporal couplings are defined by the relations 
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5^ 
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(14) 
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respectively. Then Eqs. (IT21) . ([T3| . (1141) provide us with 
following expressions for the energy and pressure [l9| : 




1 d^xP = 




[(at 
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where the generalizations of the /3-function are: 
dg(pt,fJ-s) 



(16) 



/3a(5(/i)) = 



9 In fia 



(17) 



The first term in the curvy brackets in each of Eqs. ([T5]) 
and pB]) is a classical contribution coming from the clas- 
sical energy- momentum tensor ([TJ. The second term in 
these expressions corresponds to the quantum correction 
related to the conformal trace anomaly ([2]). Thus, it 
is natural that the first terms cancel in the trace of the 
energy-momentum tensor, 7^^ = e— BP, while the second 
terms give us (fTO|) because of the equality (3t + (3s — P- 

In a lattice regularization, relations (II 5| and (jl6p were 
analytically calculated long time ago [17|. Explicit ex- 
pressions of the specific /3- functions ([TTl) are given in jTH . 

Gluon energy and pressure via gluon conden- 
sates in presence of heavy quarks and antiquarks. 

Now imagine that we have inserted a very heavy quark- 
antiquark pair in the a thermal vacuum of Yang-Mills 
theory. Since the quark and the antiquark are very heavy 
they cannot be excited by the thermal fluctuations. Nev- 
ertheless, the gluon field of the pair affects the thermal 
fiuctuations of gluons in the thermal bath, and, conse- 
quently, makes a contribution to energy density and pres- 
sure around the quark sources. 

How to take into account the presence of the heavy 
quarks on thermodynamics of gluons? An obvious way is 
to evaluate the gluon condensates in Eqs. (fTS]) and (fT6|) 
with insertions of quark creation operators. 

A static quark at the position x is created with the 
help of the gauge- invariant Polyakov loop operator 



l/T 



(18) 



where the operator P implies the path ordering of the 
integration which goes along the straight line parallel 
to compactified (temperature) direction of the Euclidean 
space-time. The temperature is introduced via the imag- 
inary time formalism. The length of the compactified 



direction is l/T. The path of the integration in Eq. (fT8)) 
is closed via the periodic boundary condition imposed 
in the compactified direction. The conjugated operator 
L^(af) creates an antiquark at the spatial point x. 

A system Q = Qi ■ ■ ■ Qnq Qi ■ ■ ■ Qnq, consisting of 
Nq quarks located at spatial positions Xi and Nq anti- 
quarks at points ijk, is created by the operator 

Lq{{R}) = L(fi) . . . L{xM^) • Lt(yi) . . . LHmo) , (19) 

where {R.,} = {xi, . . . xnq , yi, • . . VNq}- Since general- 
izations to various representations is straightforward, we 
consider here the color-averaged operators (fTO|) only. 

According to Eqs. ([TS]) and ([TC|) . in order to estimate 
the effect of the heavy (anti)quark system on the thermo- 
dynamics of gluons one should evaluate the chromoelec- 
tric and chromomagnetic condensates in the presence of 
the heavy quarks. Since Eqs. (jisp and ([T6|) are linear in 
the condensates, it is convenient to consider the excess, 
respectively, in the energy and in the (integrated-over- 
volume) pressure caused by the presence of the heavy 
quarks: 



^^Pe(T) = 



d^x [eaix) - e{x)] , (20) 
d3x[Ps(x)-P(a;)] . (21) 



Notice, that both ([^0]) and (PT|) do not depend on the spa- 
tial size of the system. These quantities are finite but, 
in general, nonzero. We used the symbol V (reminiscent 
of V-P) to denote the integral over the spatial volume, 
jd'^x. In order to avoid an abuse of notations, hereafter 
we remove the (5-signs in front of excesses in thermody- 
namic variables, so that Eq = SEq, Wq = 5Wq etc. 
The contribution of the quark system Q into both 



gluon condensates O = 
(0>Q,T 



/ d^x E"^ and O 



Jd^xB^ 



(22) 



{Lq)t 

where the product Lg of the Polyakov loops is defined in 
Eq. dUl). The right hand side of Eq. ([22]) is independent 
on the imaginary time. 

Substituting (IT5|) and ^ into Eqs. ^ and dH]) we 
get the excesses in the pressure and the energy of the 
system due to the presence of the heavy (anti) quarks 



EqIT) 




(23) 



(24) 
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Notice that the expectation value {{0))q^t in (ITSt . 
gets automatically replaced by {0)q^t in ((23)) . (l24l) . 

Equations ([^ and ([Ml) can already be used for a 
numerical estimation of the contribution of the heavy 
(anti)quark systems into the gluon thermodynamics. 
However, there is an analytical way to proceed further. 

Generalization of Michael- Rothe sum rules. 

Consider a quark Q and an antiquark, Q, separated by 
the distance R. The color-averaged QQ potential Vqq is 
given by the expectation value [20] : 



(25) 



Its important to notice that the thermodynamics of 
heavy quarks is often discussed with the help of the renor- 
malized Polyakov loops U°'^{R) that are multiplied by a 
perimeter/length-dependent factor Z: 

L{R) ^ L-"(i?) = [Z{gr^^ ■ L{R) . (26) 

However, at this stage the renormalization is irrelevant 
for the determination of both the energy p3|) and the 
pressure ([M)) . because Eq. (|22l) ~ used both in and 
in (|24p - is invariant under the renormalization shift (1261) . 

The contribution of the finite-sized heavy-quark sys- 
tem into the thermodynamics is not an extensive quan- 
tity. Thus, the heavy quark potential (|25|) is a function 
of four dimensionfull parameters: the temperature T, the 
distance between the quark and antiquark R, and the 
spatial and temporal renormalization scales and fit- In 
the presence of the heavy quarks we rewrite the dimen- 
sional identities ()13|) and the scale-independence require- 
ments ini) using the quantity ln(L(0)Li'(^)) ee -V, 
instead of InZ. The first pair of equations is: 



d 



\dlnR ain/i^/ T 
/ d d \ VQQ{R,T;^i„nt) 

\d\nT dlniij T 
while the second pair gives us 

d Von 



0, 



0, 



dT T 



/'■I 


la 







R- 



dVr 



QQ 



dR 



/'■■I 


1(51 







B\x)~E'{x) 



E^{x) + B^{x) 



B^{x)-E^ix) 



E^{x)+B^{x) 



QQ-.T 



(27) 
(28) 

(29) 



QQ,T 



(30) 



9 

Equations (|29|) and ([30| are new finite-temperature 
sum rules that relate the QQ potential to the excesses 
in expectations values of the chromoelectric and chromo- 
magnetic condensates, that arise due to the presence of 
the external quark sources. 



These sum rules is a generalization of already known 
low-energy relations that were derived in T = lattice 
gauge theory in Refs. For example, subtracting 

Eq. ([291) from Eq. ([30]) and using ([18|) one gets: 

Z^^^-*^^ U'xCix)) (31) 



1 



"In--' 



grp j ^QQ - 
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' QQ,T 

This relation is nothing but a natural finite-temperature 
extension of a well-known "action sum rule" , 

^„,H-fl^.(?M/A£W) . (32) 

^ 9 J ' QQ{T=0) 

Moreover, a zero-temperature limit of (12^ gives us a con- 
tinuum version of the so-called "energy sum rule" : 



Vnn 



d^a;. 
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B\x)~E^{x) 



E^{x)+B^{x) 



QQ(T=0) 



QQ(T=0)_ 



(33) 



Let us mention that we have three pairs of equations 
that look similar but have very different meanings: 

i) The pair (|15p and (jl6p is an equality, that relates 
the energy and pressure of the gluonic fields to the 
expectation values of the gluon condensates. This 
pair is a continuum version of the corresponding 
lattice formulae derived in Ref. [17| . 



ii) The pair (|23p and ([M)) is a natural definition of the 
contribution of heavy quarks to the energy ([T5|) and 
pressure ([TB| of the gluons. 

iii) Finally, the pair ([29| and ([30]) describes the new 
finite-temperature sum rules. 

New sum rules and exact relations for the 
heavy-quark thermodynamics. A comparison of the 
exact relations ([23|), ([Mj) with new sum rules ([50)) 
shows that the excess in the energy ([20)) and the (volume- 
integrated) excess in the pressure ([2T)) due to the presence 
of the heavy quark- antiquark pair are, respectively: 



Eqq{R,T)^Vc 
"PQQiR^T) 



QQ~ ^ ~gx~ ^ ^ 
RdVQQ 



d Vqq 
dT T 



, (34) 
(35) 



These relations are exact. 

From Eqs. ()34l) and ([55)) we get a Maxwell-type relation 
between the excesses in the energy and pressure: 

where the subscripts mean that the derivatives are taken 
at the fixed temperature T and the fixed interquark dis- 
tance i?, respectively. Equation ()36|) is an exact equation 
of state for heavy quarks. 



R 
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A comparison of the energy excess (p4|) with the cor- 
responding thermodynamical relation [the first formula 
in pT|) ] would lead us to the conclusion that the heavy 
quark potential Vqq may be interpreted - up to a renor- 
malization constant (pS)) - as an excess in the Helmholtz 
free energy Fqq due to the presence of the heavy quarks. 
However, the excess in the pressure ([35]) has nothing to 
do with its analogue in the thermodynamical limit [the 
second formula in making the free energy interpre- 

tation of the heavy quark potential Vqq obscure. 

Nevertheless, one can still interpret Vqq as an excess 
in the Helmholtz free energy of gluons due to the pres- 
ence of the heavy (anti)quarks. Indeed, besides the tem- 
perature T, the heavy quark thermodynamics has the 
additional external variable R but lacks the usual vol- 
ume variable V 21|. In this case the usual fundamental 
thermodynamic relation, dE — TdS — PdV, should be 
written as follows 



(37) 



where Sqq is the excess in the entropy and Xqq is the 
generalized force associated with the distance R between 
the quark and antiquark. From Eq. p4p one gets: 



Snn — 



QQ 



V dT 



R 



(38) 
(39) 



so that the generalized force Xqq is related to the (inte- 
grated) pressure (1551) . Then the differential of the excess 
in the Helmholtz free energy Fqq = Er 



d-pQQ - -Sqq - Xqq dR . 



(40) 



This formula is consistent with other relations: the pair 
of relations pOt . (|38)) restores the energy relation (l34t . 
while the pair (l40l) . (|39p gives us back the pressure ex- 
cess ([55)) . 

The excess in the free energy Fqq can be easily ob- 
tained by an integration of Eq. (l40l) : 



(41) 



Fqq{R,T)^Vqq{R,T). 



A free integration constant - which inevitably appears 
in any integration - is to be interpreted as the renor- 
malization quark-antiquark potential (1261) . We have just 
shown that Eq. (|4T|) - despite its standard appearance 
- has somewhat nonstandard interpretation because it 
involves the generalized fundamental thermodynamic re- 
lation (|37l) that includes the generalized force . 

The excess in the pressure due to the presence of the 
heavy quark pairs is always negative because the 
Vqq potential is a concave function of the distance (23j . 
Thus, the presence of the quarkonium states should de- 
crease the pressure in the quark-gluon plasma. A gas 



of weakly-interacting heavy-quark bound states with the 
density p should produce the total pressure deficit 



piR,T)RdVQQ{R,T) 

AFqq{R,1)- ^ 



<0, (42) 



(we assumed that all QQ-states have the same size R). 

It is interesting to notice that even at zero tempera- 
ture the integrated pressure (|35)) is not a positively de- 
fined quantity. This fact does not contradict the com- 
mon wisdom. For example, a similar feature character- 
izes the well-known Casimir effect: the vacuum between 
two perfectly conducting metallic plates has a negative 
pressure, negative free energy and negative internal en- 
ergy [2J|. Less exotic example of a similar system at 
nonzero temperature is a solution of ethanol in water, 
that demonstrates in normal conditions an effect of vol- 
ume contraction (i.e., negative excess of volume). The 
maximal volume contraction is reached at approximately 
40% fraction of ethanol, thus mimicking, in very loose 
terms, a noticeable (excess of) negative pressure. 

It was noticed recently [i^ that the contribution of 
the heavy quarks to the specific heat of the gluon plasma 
becomes negative at high enough temperatures. It was 
concluded that the heavy quark potential does not be- 
have as a true thermodynamic free energy [2^ indicating 
that one should call Vqq as "heavy quark interaction en- 
ergy" . This conclusion is based on the numerical data 
for the renormalized heavy quark potential obtained in 
SU{S) lattice gauge theory [J, [26j and supported by a 
perturbative calculation '27]. We certainly agree with 
the fact [25| that one should not consider the quark con- 
tribution as the thermodynamic quantity on its own, so 
that the identification (PT|) - which we still claim to be 
correct taking into account its derivation in our paper - 
should be understood as "the excess in the free due to 
heavy quarks is equal to the (renormalized) heavy-quark 
potential" (following, e.g., [J]). Indeed, there is no ab 
initio restriction on the sign of the excess in any thermo- 
dynamic quantity due to presence of external particles. 
For example, it is well-known that adding salt to water 
reduces the specific heat of the water, exactly as heavy 
quarks should do with the plasma at very high tempera- 
tures. Nevertheless, we do not attribute a negative spe- 
cific heat to the salt thus avoiding a thermodynamical 
puzzle. 

Multiquarks and a single quark. A generaliza- 
tion of our results to multiquark systems is quite sim- 
ple. Instead of the QQ potential Vqq one should use its 
multiquark counterpart Vq defined via the multiquark 
loop (IT51) as follows: 

e-^««^>'^)/^ = (is(m)). (43) 
In addition, the derivatives with respect to the QQ dis- 
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tance should be replaced by the following scale derivative: 

• (44) 



OR 



x=i 



According to Eqs. (gS]) and (gH), a single quark Q 
cannot affect the pressure of the system even in the de- 
confinement phase. Then we get the following nonobvi- 
ous relation [an analogue of Eq. ([30|) for Q = Q]: 



9 



9 ' J ' Q.T 

which leads to an analogue of ([^^ for a single quark 

Q-T 



(45) 



Vq-T- 



dT 



(46) 



This is yet another new sum rule, that reduces to the 
known relation p2p in the limit T — > (to this end one 
should notice that dVq/dR = 0). 

A few phenomenological examples. The excesses 
in free energy (PT|) , entropy (pS)) and internal energy ([M)) 
were thoroughly studied in lattice simulations of S'[/(3) 
Yang-Mills theory 0, H^- Here we discuss briefly the 
new quantity, the (volume-integrated) excess in the pres- 
sure (|35p caused by the presence of the heavy quark- 
antiquark pair. 

In the deconfinement phase the data for the heavy- 
quark potential can be described by the relation 28] : 



e(T)T 
(RTyiT) 



where e(T) and d{T) are dimensionless temperature- 
dependent parameters, and /x(T) is the screening mass. 
Using Eq. ([35l) we find that in the deconfinement phase 
the absolute value of the QQ-induced pressure is a de- 
creasing function of the distance R: 



d + n{T)R 



Vqq{R,T)<0. (48) 



At large separations the pressure excess vanishes due to 
the screening. 

In the confinement phase the finite-T potential can 
be described at large distances (i?T 1) by the for- 
mula [111: 



Vqq =Vo+ <y{T)R + CT \n{2RT) , [T < T,] , (49) 
where Vq and C are certain constants. One gets 

CT + aR 



At zero temperature both the short and long distances 
can be described by the Cornell potential 



VQQ{R,T = 0) = --+aR 



(51) 



WQQiR,T<T, 



(50) 



with a ~ 0.471 and cr ~ 0.192 GeV^ (one can take the 
parameters, e.g., from the heavy quarkonia fits ^2^). The 
integrated pressure, 

Wqq{R,T = 0) = -^-^, (52) 

takes its maximum = -2^0^/3 ~ -200 MeV at 

the distance = a/oTct ^ 0.31 fm. 

In the confinement phase the integrated excess in the 
gluon pressure, Eqs. ([50|) and (|52|) . decreases linearly with 
the increase in the QQ separation R at large enough R. 
Physically, this effect emerges due to fluctuations of the 
QCD string that is spanned between the quark and an- 
tiquark. In turn, the fluctuations of the string cause the 
backreaction of the gluons in a close proximity to the 
pair. In fact, it is well known that the fluctuations of the 
QCD string rise with the distance between its ends i?, 
leading to a gradual widening of the string ^] and caus- 
ing a modiflcation the string potential via the Casimir ef- 
fect 31]. The large negative value of the gluonic pressure 
([50|) and (l52|) looks as a natural phenomenon because the 
standard Casimir effect is often associated with a nega- 
tive pressure (23 |. 

Summarizing, we have derived flnite-temperature 
sum rules for the internal energy ([29l) and pressure (|30|) 
of heavy quarks systems. The Michael-Rothe sum rules 
for action ([32]) and energy (|33p are recovered automati- 
cally in the zero-temperature limit. It turns out that the 
new sum rules provide us with a known relation p4p for 
the internal energy E of the quarks ((20|) . 

At the same time, these sum rules lead to the new 
expression ([35]) for the spatial integral pO)) of the ex- 
cess in the gluonic pressure, APg. The pressure excess 
in a heavy quarkonia is always negative similarly to the 
Casimir pressure in certain systems. The excess is caused 
by interactions between the heavy quarks because an iso- 
lated (anti)quark does not affect the pressure at all. 

The variation in internal energy (j37p and in free en- 
ergy ([40]) can be expressed via the entropy Sqq, Eq. ((38|l . 
and the new generalized force, Xqq (x Vg, Eq. (|39|) . 
The force Xqq comes in conjunction with the spatial size 
of the heavy quark system R, similarly to the entropy- 
temperature pair. 

We derived the exact equation of state (|36p that re- 
lates the excess in the gluonic energy to the excess in 
the gluonic pressure around the heavy-quark systems. 
The generalization to multiquark systems is straightfor- 
ward ESI), da. 



so that at large distances the pressure is a linear function 
of the separation between the quark and antiquark. 



The author is grateful to Stam Nicolis for useful dis- 
cussions. 
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